Quantum transport on small-world networks: A continuous-time quantum walk approach 
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We consider tiie quantum meciianical transport of (coherent) excitons on small-world networks (SWN). The 
SWN are build from a one-dimensional ring of A'^ nodes by randomly introducing B additional bonds between 
them. The exciton dynamics is modeled by continuous-time quantum walks and we evaluate numerically the 
ensemble averaged transition probability to reach any node of the network from the initially excited one. For 
sufficiently large B we find that the quantum mechanical transport through the SWN is, first, very fast, given 
that the limiting value of the transition probability is reached very quickly; second, that the transport does not 
lead to equipartition, given that on average the exciton is most likely to be found at the initial node. 
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I. INTRODUCTION 

Many systems encountered in nature cannot be described 
by simple lattice models. In general such systems are char- 
acterized by graphs whose bonds connect sites with a wide 
distribution of mutual distances. Examples can be found in 
various fields, ranging from physics or biology to social stud- 
ies or computer science; see |1, 2, 3] and references therein. 
More specifically, some of these systems can be described by 
small-world networks (SWN), which have large clustering co- 
efficients but short characteristic path lengths [2]. The statisti- 
cal properties of SWN have been studied to a great extent and 
are now well understood. 

A large variety of dynamical processes on graphs are re- 
lated to the spectrum of the (discrete) Laplacian of the un- 
derlying topological network 6]. For classical diffu- 
sion on SWN, which has been modeled, for instance, by ran- 
dom walks |7, 8], it was found that the probability to be still 
or again at the initial site has a complex dependence on the 
number n of steps, i.e., at short times it decays as a power- 
law of n, whereas at longer times it has a stretched expo- 
nential dependence on n. The quantum dynamics on SWN 
has been studied mainly in the framework of the localization- 
delocalization transition I9j,ilfl], where one has also assumed 
an additional (on site) disorder. Here, the transition depends 
on the complexity of the SWN. A comparison between classi- 
cal and quantum diffusion was given in Ul,1 , where a quantum 
diffusion time (defined as the time where the participation ra- 
tio of the time-dependent wave function has dropped to a cer- 
tain value) was shown to be faster than its classical counter- 
part. However, even here little consideration has been given to 
the full set of eigenvectors of such systems, which become im- 
portant in the quantum mechanical extension of the classical 
diffusion process. 

To be specific, a quantum mechanical analog of continuous- 
time random walks (CTRW) can be defined by identifying the 
Laplacian (or connectivity matrix) A of the network with the 
Hamiltonian H. For simple lattices this corresponds, in fact, 
to a nearest neighbor hopping model lfT2l [Tsl [Til ITsl [T6ll . 
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The transformation replaces the classical diffusion process 
by a quantal propagation of the excitation through the net- 
work. Due to its formal similarity to CTRW, the procedure 
was dubbed continuous-time quantum walk (CTQW). In fact, 
it is known in other branches of physics under different names, 
such as the tight-binding model in solid-state physics llTH or 
the Huckel/LCMO model in physical chemistry jist]. CTQW 
are also closely related to so-called quantum graphs (QG), see, 
for instance, lfl9il20ll2ll l221. whose connectivity matrix is de- 
fined in a similar way. However, QG explicitly consider the 
bond between two nodes in the sense that bonds may be di- 
rected and are given a varying length. Thus, CTQW are, to 
some extent, a simplified version of QG. Quite recently, Smi- 
lansky discussed the connections between discrete Laplacians 
(equivalently, between the connectivity matrices) on discrete 
QG and periodic orbits ll23ll . There is certainly a large math- 
ematical backbone on which to establish further connections, 
see, for instance, |24]. 



n. QUANTUM WALKS ON NETWORKS 

Here, we consider transport processes (CTQW and CTRW) 
on networks, which allows us to study the two extreme cases 
of transport processes on such structures, namely, purely 
quantum mechanical (CTQW) and purely classical processes 
(CTRW). Networks are a collection of N connected nodes. 
The periodicity of regular networks can be destroyed by ran- 
domly including B additional bonds into the network. In such 
a way one creates "shortcuts" and a walker can find shorter 
paths between pairs of sites than on the regular network. In 
the following we create the SWN by randomly adding bonds 
to a regular one dimensional ring, see Fig. [1] However, we 
forbid self-connections, i.e., bonds connecting one node with 
itself. 

We denote by \i) a state associated with a localized ex- 
citation at node j and take the set to be orthonormal. 
For CTRW on undirected and unweighted networks the trans- 
fer matrix is given by the (discrete) Laplacian A of the net- 
work, by which we assume equal transition rates 7 = 1 
between all nodes. The matrix A has as non-diagonal ele- 
ments Ak^j the values —1 if nodes k and j of the network 
are connected by a bond and otherwise. The diagonal el- 
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FIG. 1: Sketch of a SWN of size N 
additional bonds. 



= 16 containing B = 11 



where the index r specifies the rth reahzation of the quantity 
in question. In so doing we obtain statistical results which 
allow for a comparison with the classical ones. In particular, 
we will consider the realization-averaged transition probabili- 
ties {Trkj{t))R, the averaged probabilities {Tf{t))ji, their lower 
bound {a{t))R, and their classical analog Further- 
more, we also calculate the long time average (LTA) of each 
quantity: 



lini — 



dt ■■■ 



(5) 



ements Aj j of A equal the number of bonds fj which exit 
from node j. Quantum mechanically, the states |j) span the 
whole accessible Hilbert space; the time evolution of an ex- 
citation initially placed at node |j) is determined by the sys- 
tems' Hamiltonian H = A and reads exp(— iHt)|j), where 
we set /i = 1. The classical and quantum mechanical transi- 
tion probabilities to go from the state \i) at time to the state 
\k) in time t are given by Pk,j{t) = {k\ exp{— At)\j) and by 
TTkjit) = = I exp(-im) I respectively. By 

fixing the coupling strength between two nodes j±i | = 1, 
the time unit [fi/Hj,j±i] for the transfer between two nodes is 
set to unity. 

From the eigenvalues En of the Hamiltonian H (or Lapla- 
cian A) follows the density of states (DOS or spectral density) 
of the given system of size N, 
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The DOS contains the essential information about the system 
and shows distinct features which depend on the network's 
topology. These features also carry over to dynamical proper- 
ties, which in some cases depend only on the En- For exam- 
ple, the average classical probability to be still or again at the 
initially excited node. 
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N 
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depends solely on the En of A, but not on the eigenstates |$„) 
lUt]. In the quantum case, we find a lower bound to 7f(t) = 
Yl^j^i (0' which also depends only on the En ifisilzsll . 



7f(i) > \a{tt 



N 

E 

n=l 



(3) 



where a(t) = 4? 7^,-1 cui At). We hasten to note that the 

— — J 1 I L 

lower bound is exact for regular networks Hl5Lllq1 . The quan- 
tity |a(i) p given in Eq. (O has also been derived in a different 
context as being the form factor of QG lfl9ll . 

III. CTQW ON SWN 

We will analyze the general behavior of CTQW on SWN 
by averaging over distinct realizations R 

{■■■)R.^\Y.[--]r, (4) 



For the numerical evaluation we make use of the standard 
software package MATLAB. Specifically, we focus on SWN 
of size N = 100 with S 1, 2, 5, and 100 additional bonds; 
the ensemble average is, in general, performed over R = 500 
realizations, which guarantees a sufficiently large number of 
samples under manageable computing times. 



A. Random matrix theory 



Before going into the details of our analysis, we like to 
point to the differences and similarities of SWN with other 
approaches to study quantum transport processes. Classical 
transport over SWN differs from that over other systems, such 
as regular lattices or fractal networks, it that the transport gets 
to be faster: While the probabiUty to return to the origin de- 
cays as for regular networks, it decays as a stretched 
exponential for SWN H H, vide infra Fig. 13 a). While the 
classical dynamics over SWN is by now well-understood, lit- 
tle is known about the quantum dynamics on such networks. 

In general, several dynamical properties of networks de- 
pend only on the DOS of the system's Hamiltonian ll26ll. We 
choose the additional bonds of our SWN randomly, thus the 
corresponding Hamiltonian will have entries at random posi- 
tions in the matrix. This has to be distinguished (to some ex- 
tent) from random matrix theory (RMT) |27]. However, there 
are also similarities between RMT and SWN. The DOS of 
SWN have been compared to RMT in ll28ll . where it was found 
that the level spacing l\E = {En+i -En) of the DOS of SWN 
can be fitted by the so-called Brody distribution, which in- 
terpolates between Poissonian and Wigner-Dyson level spac- 
ings statistics, see f^S*] for details. The SWN considered in 
Ref. |28] is a Watts -Strogatz network, obtained by randomly 
permuting the bonds of a regular one-dimensional network. 
The eigenvalue statistics of random networks have been stud- 
ied in Ref. | |29l] and in the works referenced therein; the quan- 
tum dynamics on regular disordered networks has been con- 
sidered in Ii30l1 . 
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FIG. 2: (Color online) DOS p{E) (a)-(d) and level spacing distribution P{AE) (e)-(f) of SWN with N = 100 nodes and B = 1 [(a),(e)], 
2 [(b), (f)], 5 [(c), (g)], and 100 [(d), (h)] additional bonds. The lower panels (e)-(g) show also the Poissoirian (dashed line) and Wigner-Dyson 
(dashed-dotted line) statistics, panel (h) shows fits of the tails of P{AE) with different exponentials. 



Now, the DOS of a SWN differs from that of networks 
whose sites have been totally randomly connected; the DOS 
of the latter networks follow Wigner's semicircle law. Figure 
|2] shows for SWN with N = 100 nodes and B = I, 2, 5, 
and 100 additional bonds histograms of the (average) DOS 
p{E) and of the level spacing distribution P{AE), where 
AE is normalized in such a way that the average AE = 1. 
While for small B the DOS barely changes, the level spacing 
distribution shows more drastic changes, see Fig. |2] (a)-(c). 
The appearance of large isolated eigenvalues results in a non- 
vanishing P{AE) for large AE. In Figs. |2](e)-(h) [plots of 
P{AE)'\ we also show the Poissonian [exp(— Ai?), dashed 
line] and Wigner-Dyson {2r(3/2)2A£;exp[-r(3/2)2A£;2], 
dashed-dotted line} statistics. While P{AE) roughly follows 
the Poissonian statistics for B = 1 [Fig. |2] (e)], this is not 
the case when increasing B. Especially the tail of the distri- 
bution P{AE) is better fitted by the Wigner-Dyson statistics 
[Figs.|2](f) and (g)]. However, when increasing B to the or- 
der of N [Fig. |2](h)], the tail of P{AE) neither decays as 
exp(— Ai?) (dashed line) nor as exp(— Ai?-^) (dashed-dotted 



line), but rather as exp(— Ai?^), with /i w 1.2 (solid line). 
Thus, the complexity of the DOS of SWN (compared, e.g., to 
the semicircle law) leads to dynamical properties of the SWN 
not all of which can be captured by RMT. 



B. Transition probabilities 

The ensemble average of the transition probabilities 
{'^kj{t))B. allows a first glimpse on the behavior of CTQW 
on SWN. Figure [3] shows {-Kk]{t))B. for several SWN with 

= 100 nodes and different B. Note that due to the ensem- 
ble average we can choose the initial node j freely, and we 
thus take j = 50. In the absence of any additional bond, the 
excitations travel along the ring and interfere in a very regular 
manner, producing discrete quantum carpets 1 14J. Typical for 
these carpets is that they show, depending on TV, full or partial 
revivals at specific times llT4ll . 
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FIG. 3: Time dependence of the averaged transition probabilities {nkj for SWN of size A*" — 100 with (a) B = 1, (b) B = 2, (c) B = 5, 
and (d) B = 100. The initial node is j — 50 and the number of realizations is i? = 500. 



For SWN the situation is quite different. Already a few 
additional bonds obliterate the quantum carpets; the patterns 
fade away. By adding more bonds, only the initial node re- 
tains a significant value for (tt^j (t))/f at all times t. Fur- 
thermore, already for SWN with as little as S = 5 the pat- 
tern of {iTjj {t))ji becomes quite regular after a short time, see 
Fig. [S^c). This almost regular shape is reached very quickly 
when B gets to be comparable to [Fig. [3jd)]. We note, 
however, that particular realizations may still show (depend- 
ing on their actual additional bonds) strong interference pat- 
terns. These features are washed out by the ensemble average, 
so that only the dependence on the initial node stands out. 
We will return to the discussion of the transition probabilities 
(^fc, (t))flin Sec. HUD] 

For the ring the LTA can be calculated analytically. De- 
pending on whether TV is even or odd, the LTA are slightly 
different ifTill . For even N (superscript there are two max- 
ima at fc = j and at fc = j + N/2, both having the value 

Xlj = limr^cx^ ^ Jo dt nk,j{t) = {2N~2)/N^; this is due 
to the fact that the number of nodes from j to j + N/2 is the 
same in both directions, which leads to constructive interfer- 
ence. On the other hand, for odd N (superscript °) there is 
only one maximum at fc = j, xl j — — 1) /N'^. 



Figure |4] shows {xk,])R for SWN of size N = 100 with 
B = 1, 2, 5, and 100. For B = 1 and fixed j, the two 
peaks of the regular network turn into a main peak and into a 
much weaker side peak at fc = j + N/2. This structure is still 
(barely) visible for B — 2. Akeady for B = 5 the side peak 
has practically vanished; see Fig.|4fc). While for 5 = 1,2 and 
5 also structure around the main peak is visible, for B = 100, 
the {xk,j) R are sharply peaked at fc = j. We stress that this 
should not be confused with the Anderson localization, since 
there is a non-vanishing probability to go from node j to all 
other nodes k ^ j. The sharp peak of {iTjj (t)) ^ at the initial 
node j is only the result of ensemble averaging. 



C. Return probabilities 

Since CTQW on SWN always carry the information of their 
initial node j, the averaged probabilities to return to j are a 
good measure to quantify the efficiency of the transport on 
such networks fzsl]. 

Figure |5] shows in double-logarithmic scales the ensem- 
ble averages {p{t))u, (w{t))j^, and (a{t))]i for SWN with 
N = 100 nodes and 5 = 1,2, 5, and 100. For classical trans- 
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FIG. 4: Long time average {Xk,j)R for SWN of size N = 100 with 
(a) B = 1, (b) B = 2, (c) B 5, and (d) B = 100. The number of 
realizations is J? = 500. Dark regions denote large values of (xfc , j ) fl 
and bright regions low values of (xfc,j)fl- 



port [Fig.|5fa)] the initial decay of {p{t))ii occurs faster for 
larger B. The decay at intermediate times follows a power-law 
(t"^/^) for the ring (as is clear from the linear behavior in the 
scales of the figure) and changes to a stretched exponential- 
type when B is large 101 • Thus, a classical excitation will 
quickly explore the whole SWN, so that it will occupy each 
site with equal probability of already after a relatively 
short time, see the final plateau in Fig.|5ja). 

Quantum mechanically, however, the situation is more 
complex. Let us start with the ensemble average (7f(t))i^, 
shown in Fig. |3b). For a ring of TV nodes and for times 
smaller than roughly N/2 (7f(t))/j displays a quasiperiodic 
pattern (black curve), the maxima of which decay as t^^. At 
longer times interference sets in and leads to an iiTegular be- 
havior at times larger than N/2 (2^. Now, for SWN, as long 
as B is considerably less than N, the periodic pattern still re- 
mains visible; in Fig.|5jb) one can follow how an increase in 
B (red, green, and blue curves) is smoothing out the curves, 
so that both the heights of the first maxima and the depths of 
the minima decrease. At longer times the SWN patterns are 
flattened out and {n{t)) r tends towards a limiting value. With 
increasing B this asymptotic domain is reached more quickly. 
To emphasize this point we display in Fig. |6]in an enlarged 
scale the data of Fig.|5lb) in the time interval [1, 100]. Clearly, 
for larger B the crossover from the quasiperiodic behavior at 
short times to a smoothed out pattern at longer times is shifted 
to smaller t. 

In Fig. |5lc) we plot the lower bound of Tf{t), namely 
{\a{t)\'^)ii averaged over the realizations. We note that the 
overall behavior of Figs. |3b) and|5lc) is quite similar How- 
ever, the limiting values at long times differ. For the LTA of 
(7f(i))fl we have (see also Eq. (17) of Ref. 1|31(1 ) 

{x)r = ( lim i / dtn{t)) 

^ ^(^n^r - En',r)\{j\'i>n,r){j\'i>n',r)f^) 

r,j,n,n' 




FIG. 5: (Color online) Time dependence of the averaged probabilities 
(a) (b) (7f(t))fl, and (c) (|Q(t)P)jj for SWN of size iV = 

100 with B = 1, 2, 5, and 100. The number of realizations is R — 
500. 

where 5{En,r — En'^,.) = 1 for En^r = En' .r and 6{En,r — 
En',r) ~ otherwise. For {\a{t)\'^)fj the long-time values 
for different B collapse to one value. In fact, the LTA of 

(|a(t)|2)^obeys 

(7) 

as can be immediately inferred from Eq. (O. Thus this quan- 
tity is only a function of the eigenvalues E„.r and does not 
depend on the eigenstates |$„ ,■)• In order to quantify the dif- 
ferences between Eqs. (|6|l and (|7J for SWN, we will assume 
that all the eigenvalues are nondegenerate (this assumption is, 
of course, not valid for the ring, see below). In Eq. ^ the 
triple sum adds then to RN, so that the rhs equals 1/A^. On 
the other hand, Eq. (|6]l leads to 

This expression depends on the eigenstates; in fact the rhs of 
Eq. ([8]) is the ensemble average of the averaged participation 
ratio of the eigenstates |$n,r)- Equation ^ is well known 
in the theory of quantum localization, see, e.g.. Sec. V. A. in 

iH. 
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FIG. 6: (Color online) Zoom into Fig. |5jb) for short times t — 



Now, Fig. Q shows the behavior of (x)j?, according to 
Eq. ©, for a SWN with N = 100, 500, and 1000 nodes 
as a function of B/N (we restrict ourselves to even N, the 
case of odd N is similar). Increasing B results in an increase 
of {x)r, starting from the corresponding value for the ring 
(B = 0, only one realization, and N even) 



2N ■ 



N ^ 

j 



iV2 



(9) 



where Xjj = (2A^ - 2) /N"^. From Eq. Q we obtain a 1/iV de- 
pendence for the LTA of {\a{t)\'^) ji, which by rescaling with 
(Xiing)iJ ~ 1/-^ would result in a constant value for large 
N. However, rescaling (x) r with (Xring) R shows an increase 
with N of (x)/j/(Xring)fl: which is less than linear, thus, (x) r. 
depends on A'' as l/iV", with ly e [1,2]. Additionally, for 
larger TV (see N = 500 and 1000), (x) r has a maximum 
value at B/N « 0.14, which is not present for smaller N 
(see N = 100), meaning that for this ratio of B/N the trans- 
port from the initial node to all others is least probable, a fact 
which remains unclear. A detailed study of the N dependence 
will be given elsewhere. When increasing B to the order of 
N, (x)fl' saturates to a plateau which increases monotonically 
with N. Thus, an increase in the number of nodes leads to a 
less probable transport from the initial node to all others. 
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FIG. 7: (Color online) The LTA of {Tf{t))R, {x)r, for SWN with 
N = 100, 500, and 1000 nodes as a function oiB/N. 

We further note that with increasing B the structures of 
{\a{t)\'^)ii and {Tf{t))]i differ even at short times, while for 
the ring the relation W{t) = holds exactly. 



In Ref. UM we showed that {p{t))R and {Tf{t))R (or 
{\a{t)\'^) n) can be regarded as measures for the efficiency of 
the excitonic transport. When increasing B, the initial quan- 
tum transport through the SWN takes place - on average - dur- 
ing a very short time scale (see Fig. O compared to the ring, 
where an excitation takes about t = N/ 2 to travel around the 
ring Hill . Additionally and in contrast to the classical case, 
where the limiting value is always given by the equipartition 
value 1/iV, for CTQW the limiting probability to be still or 
again at the initial node increases with B. Thus, an exciton 
is (on average) more likely to be found at the initial node, a 
feature which is not captured by the lower bound 
Therefore, [as, for instance, shown in Fig. |3c)] 

does not capture fine details of the transport, which the full 
expression (7f(t))fl does. 



D. Participation ratio of eigenstates 



For the ring the eigenstates are Bloch states, 



1 



N 



(10) 



from which |(fc|«'„)| = 1/^^^ follows for all |$„). By 
naively inserting this result into Eq. ^ one obtains (x)_r = 
1/A^, which differs from the exact result, Eq. (|9]l, by a factor 
of 2. The reason for this difference is that for a ring most of 
the eigenvalues are doubly degenerate. For SWN, on the other 
hand, most eigenvalues are non-degenerate. The fact that, as 
is evident from Fig.|2l {xjR for SWN increases with increas- 

I |4 

ing B points towards a change of the | (A; | $„ ) | from the value 
\/N'^. In order to quantify the difference to the ring case we 
plot in Fig.[8]the average distribution of eigenstates. 



(11) 



for SWN with N = 100 with B = 1, 2, 5, and 100. From 
Fig. [S] we remark that the l^n.j) r increase with increas- 
ing B. Additionally, the fluctuations between different val- 
ues of (S„.j)fl become larger, too. This results in a sub- 
stantial increase of (x) i? for larger B. We stress the par- 
ticular role played by the eigenstate |$o) = iV~^/^ 
which corresponds to the eigenvalue Eq = Q and for which 
(So.j) R = l/N^. Most of the other states contribute more to 
(x) R- In particular for SWN with large B, Fig.[8td), one finds 
large values for (S„ j)^; close to the band edges of £"„ (i.e., 
for n close to and close to N), in accordance with previous 
work; see, for instance Ref. 13311 . 

The situation may be visualized as follows: For the ring all 
eigenstates are Bloch states and hence are completely delocal- 
ized. Going over to SWN and increasing the number of addi- 
tional bonds B leads to localized states at the band edges and 
to fairly delocalized states well inside the band. The increase 
of {x)r shown in Fig. |6]is thus mainly due to the localized 
band edge states. 
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FIG. 8: The function {E,^j)r, Eq. (Ui, for SWN of size N = 100 witii (a) B = 1, (b) B = 2, (c) B = 5, and (d) B = 100. Note tiie 
different scaling of the z-axis in (d). The number of realizations is 7? = 500. 



The participation ratio also dominates the transition proba- 
bilities {-Kkjit)) R, which were presented in Fig.[3]in Sec. lIII Bl 
In general, the TTkj{t) = |(fc| exp(— iHt)|j)p averaged over 
the distinct realizations read: 



-iE„ rt 



(12) 



Under the assumption that the eigenvalues of SWN are non- 
degenerate, we obtain for the initial node j 



R 



'-'*l01$«..)l'l(j|<i>n'..)P 



(.13) 



The fluctuations for larger t [i-dependent sum in Eq. ( fTsl ll 
become suppressed due to the ensemble average. As can be 
inferred from Figs. [3a)-(c), when increasing B from B = Q 
only slightly up to B/N = 0.05, the fluctuations are already 
strongly suppressed. Larger values of B, see Fig. ^d) for 
B/N — 1, result in a very strong peak at the initial node j. 



Hence, the fluctuations at the other nodes k j become more 
and more suppressed in the ensemble average when increasing 
B. 

Now, averaging the time-independent term of Eq. ( fTsT l over 
all nodes j one recovers the LTA of (7f(t))j^, see Eqs. ^ and 
®: 



AT- 



(14) 



In the ensemble average, all nodes j can be considered 
roughly equal, thus every node j gives approximately the 
same contribution to the sum over j and we get therefore 



{X)i 



TTjj {t))u. Figure [7] shows 



that for increasing B the LTA {x)r is always larger than 
(2iV — 2) /iV^ (the corresponding value for the ring), also lead- 
ing to the almost regular shape of the transition probabilities 
{TTkj{t))R shown in Fig. [3] As noted earlier, single realiza- 
tions may still show strong interference patterns. For QG, 
Kottos and Schanz have given conditions for finding almost 
scarred eigenfunctions (states with excess density near unsta- 
ble periodic orbits of the corresponding classical chaotic sys- 
tem) ll22ll . In combination with Smilanskys work on discrete 
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QG f23*l, it might be possible in the future to obtain similar 
conditions for the networks considered here. 

We stress again that there is no Anderson localization in our 
system. Although the states are localized for large B, there is 
still a non-vanishing transition probability to go from the ini- 
tial node j to all other nodes. Thus, the additional bonds in 
the SWN do not prohibit the transport through the network 
completely, but just hinder it. Adding disorder to our system 
will essentially result in the model considered in Ref. ifToll . 
In this work, the Anderson model was augmented by addi- 
tional bonds, such that a SWN develops, which lead to the 
locaUzation-delocahzation transition. 



IV. CONCLUSION 

We modeled the quantum mechanical transport of (coher- 
ent) excitons on small-world networks by continuous-time 
quantum walks and computed the ensemble average of the 



transition probability to go from one node of the network to 
any other node. The transport through the network turns out 
to get faster with increasing the number of additional bonds. 
Distinct from the classical case, however, where the informa- 
tion of the initial node is quickly lost, quantum mechanically 
this information is preserved. During its time development 
the exciton is on average most likely to be found at the ini- 
tial node. The reason for this is to be found in the network's 
eigenstates, which are localized at the band edges, whereas 
they are quite delocalized inside the band. 
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